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Abstract 

We apply the linear delta expansion to the quantum mechanical version of the 
slow rollover transition which is an important feature of inflationary models of the 
early universe. The method, which goes beyond the Gaussian approximation, gives 
results which stay close to the exact solution for longer than previous methods. It 
provides a promising basis for extension to a full field theoretic treatment. 
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1 Introduction 

Inflationary models of the early Universe rely on the slow evolution of an inflaton fleld 
If from the initial unstable vacuum state in which {(f) = to the flnal stable vacuum in 
which {(f) = ±a, say. The effective potential Ves{(f>c) giving rise to this transition has the 
generic form of a gentle hill centred at (fc = with minima at (fc = ±a. 

The transition can be discussed at various levels of sophistication. At the most naive 
level one can think classically in terms of a ball rolling slowly down the slope of the poten- 
tial. The corresponding quantum-mechanical problem, which is the subject of the present 
paper, is the time-development of a state whose wave function is initially concentrated 
around the position of the maximum of the potential. The full treatment of the problem 
must, of course, be formulated within the framework of quantum field theory. 

The first treatment of the quantum mechanical problem was given by Guth and Pi 
who solved exactly the equation of motion for an initial Gaussian wave- function in 
an upside-down harmonic oscillator potential V = —^kx^. This was followed by a paper 
by Cooper et al. [Q, who used a Gaussian ansatz in the Dirac time-dependent variational 
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principle for the standard symmetry-breaking potential V = A(x^ — a^)^/24. The resulting 
"Hartree-Fock" solution tracks the exact solution for a short time, but departs from it 
before the time at which (x^) reaches its first maximum. Several years later Cheetham and 
Copeland went beyond the Gaussian approximation by using an ansatz which included 
a second-order Hermite polynomial. This represented an improvement on the Hartree- 
Fock approximation, but still did not reproduce the first maximum in (x^) of the exact 
wave-function. 

In the present paper we tackle this problem afresh using the linear delta expansion. 
This is a method akin to perturbation theory, but with the crucial difference that the form 
of the unperturbed Hamiltonian Hq is not fixed once and for all, but varied at each order in 
the expansion by some well-defined criterion. The role of the formal parameter S is simply 
to keep track of the order of the expansion. The method has the great advantage that 
its generalization to field theory is straightforward. We are motivated to apply it to this 
dynamical problem by its success with the static properties of the anharmonic oscillator, 
where it can be proved rigorously that it converges to the exact result when applied to the 
finite-temperature partition function|T|] and the energy levels]^]. 

The relevant Hamiltonian {h = 1) is 

1 

H = — — + \{x^ - 0^)^/24: + const. 



2 9x2 

1 92 1 



-m X + gx , (1) 



2 9x2 2 

with m2 = and g = A/24, which we split according to 

H = -\^^±lf,^x' + 6g{x'-px% (2) 

where 2gp = w? =p p^. That is, we choose as bare mass term ±|/i2x2. The sign of the 
term as well as the value of p will be determined as functions of t after the perturbative 
expansion has been carried out to a given order (at which stage 5 is set equal to 1) by the 
criterion of minimal sensitivity (PMS) 0, namely that 

Note that for either sign of the new mass term, /i has a limited range. In case (i), when 
the mass term is — we have 2gp = th? — p?. The essence of the delta expansion is that 
the extra term —gpx"^ in the interaction should compensate as far as possible the original 
term gx^., which means that p should be positive. Hence we require that p? < rn?. In 
case (ii), when the mass term is +|/^^, the same restriction will arise from the form of the 
zeroth-order solution and the initial wave-function. 
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2 Delta Expansion 



The two cases need to be treated separately. Which of them is relevant at a given value of 
t is determined by the PMS criterion. 



2.1 Case (i) 

In this case the bare Hamiltonian is 

1 ^2 1 



It is useful to scale x according to x — y/^/ji, so that 



Ho = -^l, + yA. (6) 



2 \dy^ J 

Given that the initial wave-function is a Gaussian of the form i\}(t = 0) = Aexp(— i^y^), the 
zeroth-order equation of motion i^oV^o = '^c^V'o / can be solved exactly by a wave- function 
of the same form with A and B becoming functions of t. The equations they have to satisfy 
are 

iA/fx = AB, (6) 

with solutions 

S = ^ tan(?7o - ifit) 

A = J\f (cos(rio - iiit)y^ (7) 

where riQ is determined by B{t = 0) = ^tanr^o and the normalization constant J\f by 

A{t = Q) = N (cos?7o)~^. This is precisely the solution of Guth and Pi for the upside-down 
oscillator, with m replaced by /x. 

To obtain a systematic perturbative expansion in powers of 8 it is useful to write 
■0 = </?exp(— S(i)y2). The equation for (/? is then 



Bip + 2Byip'-^ip" + S~g{y'-m'] 



(8) 



where we have scaled g and p according to gf = p^g and p = p/fi. 

A general feature of perturbative expansions for a polynomial potential of degree p 
is that (fi is also a polynomial, of degree Np in A^th order of the expansion. Thus in the 
present case, expanding (p ets (p — J2 ^^fn, the first-order part (fi is an (even) polynomial of 
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degree 4, which we write as (fi — a + by^ + cy'^. The equations of motion for the coefficients 
a, b and c are 



ia/fj, ~ Ba — h 

ib/ II — 5Bb — 6c — pgA 

id/ II — 9Bc + gA, 



(9) 



which can be solved successively in reverse order, using the solutions for A and B previously 
determined. The initial conditions ai t — are a — b — c — 0. Thus c is given by 



igAf/8 |3g^2sinh2^+^sinh4g-co|, 

(cosh^)9/2l 4 



(10) 



where 9 = iit + irjo and Cq = Sirjo + 2isin2?7o + (V^) sin 4770- 
Using this solution in the equation for b we obtain 

—ioKf rl~l ~ 3ir~~ 

b= ^^^i — + - sinh 2^) + 60 + -r 3^ tanh ^ + cosh^ ^ - cq tanh ^ 

(cosh^)^/2 1 2 2 4 1- 



where 



X % 9i 3 3i 

^0 ^ ^PVo + |Psin2r7o + — r^otanr^o - ^Cotanr/o - — cos^r^o- 



Finally, using this solution in the equation for a we obtain 

a — — — I — -5^tanh6' + 6otanh6' — On 

(cosh^)V2\ 2^ 



+ 



3i 



■l-9sech'^9 + ^ tanh §) + 9 + ^cosech^^i 



(12) 



where 



oo = -prio tan?7o + ibo tanr/o + -i 



(^-i?7osec^?7o + i tan r/o) + iTjo + -cosec^r/o 



2.2 Case (ii) 

The zeroth-order equations in this case are 



iB/ii 
iA/n 



2B' 
AB, 



2' 



(13) 



with solutions 



B — -coth(?7o + ipt) 

A = N{smh{riQ-\- in,t))~^ . 



(14) 
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As mentioned in the Introduction, the restriction on /x in this case comes from the form 
oi B{t = 0) and the form of the initial wave-function, which, in all the papers quoted, 
is taken as a minimal wave-packet appropriate to a positive mass term +|m^a;^. In the 
present formulation this means that B{t = 0) = (l/2)(m//i). But since B{t = 0) = 
(l/2)coth?7o < 1/2, we have the same restriction on /i, namely /i < m, as in Case (i). 

The first-order equations for a, b and c are identical in form to Eq. (^), but the driving 
terms A and B are now different. 

The solution for c is now 



-^HZ^(3^~-2sin2^~+isin4^~ 
(zsin^)9/2l 4 



Co 



(15) 



where Cq = —SirjQ + 2i sinh 2?7o — (V^) sinh4?7o and in this case 6 = [it — itiq. 
Using this solution in the equation for h we obtain 



, , '-p(e--sin2^) + 6o-- 
(^sin^)5/2 L 2^^^ 2 ^ " 4 



-36'cot6' + cos^ e + cocot6' 



(16) 



where 



1 , 9 , 3i , 3,2 

jP^o - smh 2?7o - -rjo coth rjQ + — Cocothr^o + - cosh rjQ. 



Finally, using this solution in the equation for a we obtain 



(i sin ^) 1/2 



-pOcotO + ibocotO + ao 



3i 



3~ 

-^cosec^6' + -cot6' — 9 — -cocosec^O 



(17) 



where 



«o = 2^^ocoth?7o 



6ocoth?7o - - 



^?7ocosech^?7o + ^coth?7o + r/o - ^cocosech^?7o 



We have checked these solutions by numerical integration using the Runge-Kutta method. 
This reveals that in Case (ii) care needs to be taken to ensure that we are on the appropriate 
branch of the square roots. At values of t where sin^ = — 1, a naive numerical evaluation 
will stay on the first sheet, thus giving rise to a discontinuity, whereas the true solution is, 
of course, continuous. 

In fact, as we shall see, the coefficient a is not needed in the calculation of (x^)2 to first 
order in 6, though it would, of course, be needed in higher order. 

3 Variational Aspect 



The expressions we have obtained all depend on the parameter p introduced in Eq. (0). 
The other essential aspect of the delta expansion is that such a parameter is determined 
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by some non-perturbative criterion, most frequently the principle of minimal sensitivity, 
Eq. (I). 

To that end we need an expression for (a;^), which, given that the wave-function is a 
(complex) Gaussian with polynomial corrections, can be written down in closed form in 
terms of the coefficients A,B,a, b, c. Thus to order 6, 



IV'I' = [\A\^ + 26Re {A*{a + by^ + cy^)} 
where a = 2ReB, so that 

1 1 + {26/\A'^\)Re{A*{a + 3b/{2a) + 15c/ (4a2))} 



if 9A f 7i 11 

(19) 



2a 1 + {26/\A^\)Re {A*{a + b/{2a) + 3c/(4a2))} 
1 

A^ a a^' 



2a 



It is an interesting feature of the structure of the perturbative equations that the wave- 
function is automatically normalized to the order we are working. That is, / (p^^i = 0. 
Thus the second equation is identical to the first, and not merely an 0(5) approximation 
to it. The expectation value we seek is obtained on scaling by /i, i.e. (x^) = 

At this stage we set 5 = 1 and apply Eq. (^. This has to be done for each time t, 
and the result is that the chosen value /i of /x now becomes a function of t, even though 
H was treated as a constant in the equations of motion. In the present case, since we are 
unable to go to very high orders in the expansion this is a more important property than 
the iV-dependence of /i. The 0(5°) calculation does not have such a stationary point. 

In Fig. 1 we show graphs of (x^)2 for various values of t. The parameters chosen are 
those used in Refs. Q and [^], namely a = 5 and A = 0.01 (which corresponds to a "large" 
dimensionless coupling contant^]). We include both cases by plotting (x^)2 as a function 
of (T/i, where a = —1 for case (i) and -|-1 for case (ii). There is a well-defined maximum 
which moves steadily to the right as t increases, crossing over from case (i) to case (ii) at 
about t = 11. From these and similar graphs we extract the value of fl{t), which is plotted 
in Fig. 2. 

Using these values of /2(t) we can then calculate (x^) 2 (/2(t), t) from Eq. ([T9|) as a function 
of t. This is plotted in Fig. 3 along with the results obtained using the "Hartree-Fock" 
method of Ref. 0, the improved variational method of Ref. 0, the exact value of (a;^)2 (t), 
obtained by Fourier transform and numerical integration 0, and finally the result of first- 
order perturbation theory. The latter corresponds to 0(5) of the delta expansion, but with 
fixed at m in case (i), and exemplifies the importance of the t-dependence of ft. 

As can be seen, the delta-expansion result tracks the exact result for longer than either 
of the other variational calculations, essentially up to the point where (x^)2 reaches its 
maximum, but then overshoots. A similar degree of accuracy in quantum field theory 
would mean that, to this order of the expansion, the infiationary period would be very well 
described, but the reheating process less so. To extend the range of the approximation to 
longer times a higher-order calculation would presumably be needed 
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4 Discussion 



Figure 3 is our main result, but it is also of interest to enquire how closely the calculated 
wave-function agrees with the exact result, since a well-known feature of variational meth- 
ods is that quite reasonable values for expectation values such as (x^) can be obtained 
with rather inaccurate wave-functions. In fact our wave-function agrees rather well with 
the true wave-function up to t ~ 6, but begins to diverge from it thereafter, even though 
still giving good values for (x^). In Fig. 4 we plot the two values of versus x for t = 6 
and 8. 

Various extensions of the present treatment are possible. Given the simplicity of the 
first-order equations resulting from the method it seems that the extension to 0(5^) should 
be relatively straightforward, certainly if the integrations are performed numerically, and 
it would be interesting to see the improvement thereby achieved. As mentioned earlier, the 
wave-function would involve an even polynomial of order 8 multiplying the lowest-order 
Gaussian. Another possibility is the use of the original delta expansion^, whereby the 

term in the potential is written as x'^^^^^^ and expanded as + Slnx^ + . . .). This 
expansion is known to converge for the energy levels of the anharmonic oscillator, and the 
0{6) calculation for the present problem should be tractable. However, the disadvantage 
of this method is that its extension to field theory beyond first order becomes extremely 
difficult. 

The most important extension is clearly to attempt to apply the methodology of the 
linear delta expansion to the quantum field theory problem. The importance of going 
beyond the Gaussian approximation has been emphasized in refs. and f^], and technically 
the linear delta expansion is essentially a modified perturbation theory, modulo the crucial 
variational aspect. 
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Figure 1: Graphs of (a;^)^ versus a/i for t—5, 9 and 13, where a = —1 for case (i) and 
a = +1 for case (ii). We have excluded the region /i < 0.05 since there are severe round-off 
problems near jj, = 0. 
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Figure 2: fj, versus t. The change-over from case (i) to case (ii) occurs between t—11 and 



10 




Figure 3: {x^)^ versus t. First-order linear delta expansion (LDE) compared with the 
exact result (Exact), the variational calculations of Ref. 0] (HF) and Ref. (CC), and 
first-order perturbation theory (PT). 
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Figure 4: Graphs of versus x for t=6 and 8. The sohd hne is the first-order LDE 
calculation and the dotted line is the exact result. 
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